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1 Introduction 



The least well quantified parameters of the Standard Model (SM) are the masses of light 
quarks and the p and r\ parameters in the Wolfenstein representation of the CKM mixing 
matrix. A non-zero value of 77 signals CP violation. The important question is whether the 
CKM ansatz explains all observed CP violation. This can be addressed by comparing the SM 
estimates of the two CP violating parameters e and e'/ e against experimental measurements. 
The focus of this talk is to evaluate the dependence of these parameters on the light quark 
masses and on the bag parameters Bk, B^ 2 , and B^ 2 . I will therefore provide a status 
report on the estimates of these quantities from lattice QCD (LQCD). 

Since this is the only lecture presenting results obtained using LQCD at this school/ workshop, 
I have been asked to give some introduction to the subject. The only way I can cover my 
charter, introduce LQCD, summarize the results, and make contact with phenomenology 
is to skip details. I shall try to overcome this shortcoming by giving adequate pointers to 
relevant literature. 



2 Lattice QCD 

LQCD calculations are a non-perturbative implementation of field theory using the Feynman 
path integral approach. The calculations proceed exactly as if the field theory was being 
solved analytically had we the ability to do the calculations. The starting point is the 
partition function in Euclidean space-time 

Z = J VA^ Vi) V~i) e~ s (1) 

where S is the QCD action 

S = jd A x (^F^-^MVO • (2) 

and M is the Dirac operator. The fermions are represented by Grassmann variables ip and 
ip. These can be integrated out exactly with the result 

Z = jvA, detM el d4x H^""). (3) 

The fermionic contribution is now contained in the highly non-local term detM, and the 
partition function is an integral over only background gauge configurations. One can write 
the action, after integration over the fermions, as S = S ga uge + S quar k s = J d A x (^F^F^) — 
J2i Ln(DetMj) where the sum is over the quark flavors distinguished by the value of the bare 
quark mass. Results for physical observables are obtained by calculating expectation values 

(O) = '-j VA, O e~ s . (4) 

where O is any given combination of operators expressed in terms of time-ordered products 
of gauge and quark fields. The quarks fields in O are, in practice, re-expressed in terms of 
quark propagators using Wick's theorem for contracting fields. In this way all dependence on 
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quarks as dynamical fields is removed. The basic building block for the fermionic quantities, 
the Feynman propagator, is given by 



S F (y,j,b;x,i,a) = (M (5) 

where M~ l is the inverse of the Dirac operator calculated on a given background field. A 
given element of this matrix {M~ 1 ) v J'? b is the amplitude for the propagation of a quark from 
site x with spin-color i, a to site-spin-color y,j, b. 

So far all of the above is standard field theory. The problem we face in QCD is how 
to actually calculate these expectation values and how to extract physical observables from 
these. I will illustrate the second part first by using as an example the mass and decay 
constant of the pion. 

Consider the 2-point correlation function, (0| J2 X ®f{%i t)Oi(0, 0) 1 0) , where the operators 
O are chosen to be the fourth component of the axial current Of = Oi = At = -^7475-^ 
as these have a large coupling to the pion. The 2-point correlation function then gives the 
amplitude for creating a state with the quantum numbers of the pion out of the vacuum at 
space-time point by the "source" operator Of, the evolution of this state to the point (x, t) 
via the QCD Hamiltonian; and finally the annihilation by the "sink" operator Of at (x,t). 
The rules of quantum mechanics tell us that Oi will create a state that is a linear combination 
of all possible eigenstates of the Hamiltonian that have the same quantum numbers as the 
pion, i.e. the pion, radial excitations of the pion, three pions in J = state, . . .. The second 
rule is that on propagating for Euclidean time t, a given eigenstate with energy E picks up 
a weight e~ Et . Thus, the 2-point function can be written in terms of a sum over all possible 
intermediate states 

(01 E £>,<*, t)o,(o)|o> = £ ( ° |O/l ^" |Ol|0> e- E '< ■ (6) 



n 



To study the properties of the pion at rest we need to isolate this state from the sum over 
n. To do this, the first simplification is to use the Fourier projection J2x as it restricts the 
sum over states to just zero-momentum states, so E n — * M n . (Note that it is sufficient to 
make the Fourier projection over either Oi or Of.) The second step to isolate the pion, i.e. 
project in the energy, consists of a combination of two strategies. One, make a clever choice 
of the operators Oi to limit the sum over states to a single state (the ideal choice is to set Oi 
equal to the quantum mechanical wave-function of the pion), and two, examine the large t 
behavior of the 2-point function where only the contribution of the lowest energy state that 
couples to Oi is significant due to the exponential damping. Then 

<0|£O,(*, t )O,(0)|0> tJoo (7) 

The right hand side is now a function of the two quantities we want since (0|Ai|7r) = M n f n . 
In this way, the mass and the decay constant are extracted from the rate of exponential 
fall-off in time and from the amplitude. 

Let me now illustrate how the left hand side is expressed in terms of the two basic 
quantities we control in the path integral - the gauge fields and the quark propagator. Using 
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Local Interpolating operators 




Extended Interpolating operators 

Figure 1: A schematic of the pion 2-point correlation function for local and non-local inter- 
polating operators. 

Wick contractions, the correlation function can be written in terms of a product of two quark 
propagators Sp, 

(0|^^ 7475 ^(x,t)V;7475^(0,0)|0) = (0\J2S F (0;x,t) lAl5 S F (x,t;0) l4l5 \0). (8) 

X X 

This correlation function is illustrated in Fig. |l|. It is important to note that one recovers 
the 2-point function corresponding to the propagation of the physical pion only after the 
functional integral over the gauge fields, as defined in Eq. |], is done. To illustrate this 
Wick contraction procedure further, consider using gauge invariant non-local operators, for 
example using O = Tp(x,t)'j4j5(Ve-'x dzwA ^)ip(y,t) where V stands for path-ordered. After 
Wick contraction the correlation function reads 

(0| X>(0; ^t)7475(Per^)5 F (f,t;y,0) 7 475(Pe^^)|0). (9) 

X 

and involves both the gauge fields and quark propagators. This correlation function would 
have the same long t behavior as shown in Eq. however, the amplitude will be different 
and consequently its relation to f n will no longer be simple. The idea of improving the 
projection of O on to the pion is to construct a suitable combination of such operators that 
approximates the pion wave-function. 

To implement such calculations of correlation functions requires the following steps. A 
way of generating the background gauge configurations and calculating the action S associ- 
ated with each; calculating the Feynman propagator on such background fields; constructing 
the desired correlation functions; doing the functional integral over the gauge fields to get 
expectation values; making fits to these expectation values, say as a function of t as in Eq. ^| 
to extract the mass and decay constant; and finally including any renormalization factors 
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needed to properly define the physical quantity. It turns out that at present the only first 
principles approach that allows us to perform these steps is LQCD. Pedagogical expose to 
LQCD can be found in |l| ^ [|, |J, and I shall only give a very brief description here. 

Lattice QCD - QCD defined on a finite space-time grid - serves two purposes. One, 
the discrete space-time lattice serves as a non-perturbative regularization scheme. At finite 
values of the lattice spacing a, which provides the ultraviolet cutoff, there are no infinities. 
Furthermore, renormalized physical quantities have a finite well behaved limit as a — > 0. 
Thus, in principle, one could do all the standard perturbative calculations using lattice 
regularization, however, these calculations are far more complicated and have no advantage 
over those done in a continuum scheme. The pre-eminent utility of transcribing QCD on the 
lattice is that LQCD can be simulated on the computer using methods analogous to those 
used in Statistical Mechanics. These simulations allow us to calculate correlation functions 
of hadronic operators and matrix elements of any operator between hadronic states in terms 
of the fundamental quark and gluon degrees of freedom following the steps discussed above. 

The only tunable input parameters in these simulations are the strong coupling constant 
and the bare masses of the quarks. Our belief is that these parameters are prescribed by 
some yet more fundamental underlying theory, however, within the context of the standard 
model they have to be fixed in terms of an equal number of experimental quantities. This 
is what is done in LQCD. Thereafter all predictions of LQCD have to match experimental 
data if QCD is the correct theory of strong interactions. 

A summary of the main points in the calculations of expectation values via simulations 
of LQCD are as follows. 

• The Yang-Mills action for gauge fields and the Dirac operator for fermions has to be 
transcribed on to the discrete space-time lattice in such a way as to preserve all the key 
properties of QCD - confinement, asymptotic freedom, chiral symmetry, topology, and a 
one-to-one relation between continuum and lattice fields. This step is the most difficult, 
and even today we do not have a really satisfactory lattice formulation that is chirally 
symmetric in the m q = limit and preserves the one-to-one relation between continuum 
and lattice fields, i.e. no doublers. In fact, the Nielson-Ninomiya theorem states that for 
a translationally invariant, local, hermitian formulation of the lattice theory one cannot 
simultaneously have chiral symmetry and no doublers ||. One important consequence 
of this theorem is that, in spite of tremendous effort, there is no viable formulation of 
chiral fermions on the lattice. For a review of the problems and attempts to solve them 

see IB 01- 

A second problem is encountered when approximating derivatives in the action by finite 
differences. As is well known this introduces discretization errors proportional to the 
lattice spacing a. These errors can be reduced by either using higher order difference 
schemes with coefficients adjusted to take into account effects of renormalization, or 
equivalently, by adding appropriate combinations of irrelevant operators to the action 
that cancel the errors order by order in a. The various approaches to improving the 
fermion and gauge actions are discussed in |9|, [10], [11]. Here I simply list the three 
most frequently used discretizations of the Dirac action - Wilson [[Tj|, Sheikholeslami- 
Wohlert (clover) [O, and staggered |TJ], which have errors of 0(a), 0(a s a) — 0(a 2 ) 
depending on the value of the coefficient of the clover term, and 0(a 2 ) respectively. 
The important point to note is that while there may not yet exist a perfect action 
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(no discretization errors) for finite a, improvement of the action is very useful but not 
necessary. Even the simplest formulation, Wilson's original gauge and fermion action 
|12| , gives the correct results in the a = limit. It is sufficient to have the ability to 



reliably extrapolate to a = to quantify and remove the discretization errors. 

• The Euclidean action S = J d 4 x{ \F (W F fa/ -TrLnM) for QCD at zero chemical potential 
is real and bounded from below. Thus e~ s in the path integral is analogous to the 
Boltzmann factor in the partition function for statistical mechanics systems, i.e. it 
can be regarded as a probability weight for generating configurations. Since S is an 
extensive quantity the configurations that dominate the functional integral are those 
that minimize the action. The "importance sampled" configurations (configurations 
with probability of occurrence given by the weight e~ s ) can be generated by setting up 
a Markov chain in exact analogy to say simulations of the Ising model. For a discussion 
of the methods used to update the configurations see (I| or the lectures by Creutz and 
Sokal in §. 

• The correlation functions are expressed as a product of quark propagators and path 
ordered product of gauge fields using Wick contractions. This part of the calculation 
is standard field theory. The only twist is that the calculation is done in Euclidean 
space-time. 

• For a given background gauge configuration, the Feynman quark propagator is a matrix 
labeled by three indices - site, spin and color. A given element of this matrix gives 
the amplitude for the propagation of a quark with some spin, color, and space-time 
point to another space-time point, spin, and color. Operationally, it is simply the 
inverse of the Dirac operator. Once space-time is made discrete and finite, the Dirac 
matrix is also finite and its inverse can be calculated numerically. The gauge fields 
live on links between the sites with the identification U^x^x + fi) = e m 9 A ^+P'/^) ; {_ e _ 
the link at site x in the fi direction is an SU(3) matrix U^x, x + fi) denoting the 
average gauge field between x and x + fi and labeled by the point x + ft/2. Also 
U^XjX — fi) = UUx — £i,x). The links and propagators can be contracted to form 
gauge invariant correlation functions as discussed above in the case of the pion. 

• On the "importance sampled" configurations, the expectation values reduce to simple 
averages of the correlation functions. The problem is that the set of background gauge 
configurations is infinite. Thus, while it is possible to calculate the correlation functions 
for specified background gauge configurations, doing the functional integral exactly is 
not feasible. It is, therefore, done numerically using monte carlo methods. 

The simplest way to understand the numerical aspects of LQCD calculations is to gain 
familiarity with the numerical treatment of any statistical mechanics system, for example 
the Ising model. The differences are: (i) the degrees of freedom in LQCD are much more 
complicated - SU(3) link matrices rather than Ising spins, and quark propagators given by 
the inverse of the Dirac operator; (ii) The action involves the highly nonlocal term Ln Det M 
which makes the update of the gauge configurations very expensive; and (iii) the correla- 
tion functions are not simple products of spin variables like the specific heat or magnetic 
susceptibility, but complicated functions of the link variables and quark propagators. 

The subtleties arising due to the fact that LQCD is a renormalizable field theory and 
not a classical statistical mechanics system come into play in the behavior of the correlation 
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functions as the lattice spacing a is varied, and in the quantum corrections that renormalize 
the input parameters (quark and gluon masses and fields) and the composite operators used 
in the study of correlation functions. At first glance it might seem that one has introduced 
an additional parameter in LQCD, the lattice spacing a, however, recall that the coupling 
a s and the cutoff a are not independent quantities but are related by the renormalization 
group 

Aqcd = -e-^/WM (/VCa))-^ + • • • , (10) 

where Aq C d is the non-perturbative scale of QCD, and P = (11 — 2n//3)/167r 2 and Pi = 
(102 — 38rij/3)/(167r 2 ) 2 are the first two, scheme independent, coefficients of the /3-function. 
In statistical mechanics systems, the lattice spacing a is a physical quantity - the intermolec- 
ular separation. In QFT it is simply the ultraviolet regulator that must eventually be taken 
to zero keeping physical quantities, like the renormalized coupling, spectrum, etc, fixed. 

The reason that lattice results are not exact is because in numerical simulations we have 
to make a number of approximations. The size of these is dictated by the computer power 
at hand. They are being improved steadily with computer technology, better numerical 
algorithms, and better theoretical understanding. To evaluate the reliability of current lattice 
results, it is important to understand the size of the various systematic errors and what is 
being done to control them. I, therefore, consider it important to discuss these next before 
moving on to results. 



3 Systematic Errors in Lattice Results 

The various sources of errors in lattice calculations are as follows. 

Statistical errors: The monte carlo method for doing the functional integral employs 
statistical sampling. The results, therefore, have statistical errors. The current understand- 
ing, based on agreement of results from ensembles generated using different algorithms and 
different initial starting configuration in the Markov process, is that the functional integral 
is dominated by a single global minimum. Also, configurations with non-trivial topology are 
properly represented in an ensemble generated using a Markov chain based on small changes 
to link variables. Another way of saying this is that the data indicate that the energy land- 
scape is simple. As a result, the statistical accuracy can be improved by simply generating 
more statistically independent configurations with current update methods. 

Finite Size errors: Using a finite space-time volume with (anti-)periodic boundary 
conditions introduces finite size effects. On sufficiently large lattices these effects can be 
analyzed in terms of interactions of the particle with its mirror images. Luscher has shown 
that in this regime these effects vanish exponentially [15| . Current estimates indicate that for 



L > 3 fermi and M^L > 6 the errors are ^ 1 %, and decrease exponentially with increasing 
L. 

Discretization errors: The discretization of the Euclidean action on a finite discrete 
lattice with spacing a leads, in general, to errors proportional to a, a™a, a 2 , .... The precise 



form of the leading term depends on the choice of the lattice action and operators ||16|| . For 
example, lattice artefacts in the fermion action modify the quark propagator M _1 at large 
p from its continuum form. Numerical data show that the coefficients of the leading term 
are large, consequently the corrections for 1/a « 2GeV are significant in many quantities, 
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10-30% \ 17\ . The reliability of lattice results, with respect to 0(a) errors, is being improved 
by a two pronged strategy. First, for a given action extrapolations to the continuum limit 
a = are performed by fitting data at a number of values of a using leading order corrections. 
Second, these extrapolations are being done for different types of actions (Wilson, Clover, 
staggered) that have significantly different discretization errors. We consider the consistency 
of the results in the a = limit as a necessary check of the reliability of the results. 

Extrapolations in Light Quark Masses: The physical u and d quark masses are 
too light to simulate on current lattices. For 1/a = 2 GeV, realistic simulations require 
L/a > 90 to avoid finite volume effects, i.e. keeping M W L > 6 where M n is the lightest 
pseudoscalar meson mass on the lattice. Current best lattice sizes are L/a = 32 for quenched 
and L/a = 24 for unquenched. Thus, to get results for quantities involving light quarks, one 
typically extrapolates in m u = from the range m s /3 — 2m s using simple polynomial fits 
based on chiral perturbation theory. For quenched simulations there are additional problems 
for m q <m s /3 as discussed below in the item on quenching errors. 

Discretization of heavy quarks: Simulations of heavy quarks (c and b) have dis- 
cretization errors of 0(ma) and 0(pa). This is because quark masses measured in lattice 
units, m c a and m^a, are of order unity for 2GeV < l/a < 5GeV. It turns out that these 
discretization errors are large even for m c . Extrapolations of lattice data from lighter masses 
to mj using HQET have also not been very reliable as the corrections are again large. The 
three most promising approaches to control these errors are non-relativistic QCD, 0(a) im- 



proved heavy Dirac, and HQET. These are discussed in |L9], ^0], |2T]. There will not be any 
discussion of heavy quark physics in this talk. 

Matching between lattice and the continuum (renormalization constants): 

Experimental data are analyzed using some continuum renormalization scheme like MS, so 
results in the lattice scheme have to be converted to this scheme. The perturbative relation 
between renormalized quantities in say MS and the lattice scheme, are in almost all cases, 
known only to 1-loop. Data show that the 0(a^) corrections can be large, ~ 10 — 50% 
depending on the quantity at hand, even after implementation of the improved perturbation 
theory technique of Lepage-Mackenzie ]TJ||. Recently, the technology to calculate these 
factors non-perturbatively has been developed and is now being exploited fl22 |. As a result, 



the reliance on perturbation theory for these matching factors will be removed. 

Operator mixing: The lattice operators that arise in the effective weak Hamiltonian 
can, in general, mix with operators of the same, higher, and lower dimensions because at 
finite a the symmetries of the lattice and continuum theories are not the same. Perturbative 
estimates of this mixing can have an even more serious problem than the uncertainties 
discussed above in the matching coefficients. In cases where there is mixing with lower 
dimensional operators, the mixing coefficients have to be known very accurately otherwise 
the power divergences overwhelm the signal as a — > 0. In cases where there is mixing, due 
to the explicit chiral symmetry breaking in Wilson like actions, with operators of the same 
dimension but with different tensor structures, the chiral behavior may again be completely 
overwhelmed by the artefacts. In both of these cases a non-perturbative calculation of the 
mixing coefficients is essential. 

Quenched approximation: The inclusion of the fermionic contribution, Ln det(M), 
in the Boltzmann factor e~ s for the generation of background gauge configurations increases 
the computational cost by a factor of 10 3 — 10 4 . The strategy, therefore, has been to initially 
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neglect this factor, and to bring all other above mentioned sources of errors under quantitative 
control. The justification is that the quenched theory retains a number of the key features of 
QCD - confinement, asymptotic freedom, and the spontaneous breaking of chiral symmetry 
- and is expected to be good to within 10 — 20% for a number of quantities. One serious 
drawback is that the quenched theory is not unitary and analysis of it shows the 

existence of unphysical singularities in the chiral limit. For example, the chiral expansions 
of pseudoscalar masses and decay constants in the quenched theory are modified in two 
ways. One, the normal chiral coefficients are different in the quenched theory, and second 
there are additional terms that are artefacts and are singular in the limit m q = |23|. 
These artefacts are expected to start becoming significant for m q < m s /3 [|25], p6| . Thus, 
in quenched simulations one of the strategies for extrapolations in the light quark masses 
is to use fits based on %PT, keep only the normal coefficients, and restrict the data to the 
range m s /3 — 2m s where the artifacts are expected to be small. In this talk I shall use this 
procedure to "define" the quenched results. 

The above mentioned systematic errors are under varying degrees of control depending on 
the quantity at hand. Of the systematics effects listed above, quenching errors are by far the 
least well quantified, and are, to first approximation, unknown. Of the remaining sources the 
two most serious are the discretization errors and the matching of renormalized operators 
between the lattice and continuum theories. An example of the latter is the connection 
between the quark mass in lattice scheme and in a perturbative scheme like MS. We shall 
discuss the status of control over these errors in more details when discussing data. 



4 Light Quark Masses from xPT 



The masses of light quarks cannot directly be measured in experiments as quarks are not 
asymptotic states. One has to extract the masses from the pattern of the observed hadron 
spectrum. Three approaches have been used to estimate these - chiral perturbation theory 
(xPT), QCD sum-rules, and lattice QCD. 

%PT relates pseudoscalar meson masses to m u , rrid, and m s . However, due to the presence 
of an overall unknown scale in the chiral Lagrangian, %PT can predict only two ratios 
amongst the three light quark masses. The current estimates are pTL GBL E9[ 





Lowest order 


Next order 


2m s /(m u + m d ) 
m u /m d 


24.2 - 25.9 
0.55 


24.4(1.5) 
0.553(43) . 



These ratios have been calculated neglecting the Kaplan-Manohar symmetry [ 30| . The subtle 
point here is that the masses /i; extracted from low energy phenomenology are related to 
the fundamental parameters mj, defined at some high scale by the underlying theory, as 
\ii = rrii + (3 * det M/rriiA x sB, where the second, correction, term is an instanton induced 
additive renormalization. For the u quark, the magnitude of this term is roughly equal to 
the xPT estimates of /i n for (3 ~ 2 pi| . Consequently, using x?T, one cannot estimate the 
size of isospin breaking from low energy phenomenology alone. At next to leading order, 



only one combination of ratios (Q' 



[mz 



— m )/ijn d — ml) as defined in 



can be 



S 



determined unambiguously from %PT. Even if one ignores the Kaplan-Manohar subtlety 
(such an approach has been discussed by Leutwyler under the assumption that the higher 
order terms are small one still needs input from sum- rules or LQCD to get absolute 
values of quark masses. 



5 Light Quark Masses from LQCD 

The most extensive and reliable results from LQCD have been obtained in the quenched 
approximation. In the last year the statistical quality of the quenched data has been im- 
proved dramatically especially by the work of the two Japanese Collaborations CP-PACS 
and JLQCD (see |32| for a recent review). Simultaneously, the lattice sizes have been pushed 



to ^ 3 fermi for the lattice spacing in the range 0.5 — 0.25 GeV -1 . In Fig. ^| we show the 
CP-PACS data obtained using Wilson fermions. To highlight the statistical improvement 
we show data at (3 = 6.0 from the next best calculation (with respect to both statistics and 
lattice size) [j33fl . 

To reliably extrapolate the lattice data to the continuum limit one needs control over 
discretization errors and over the matching relations between the lattice scheme and the 
continuum scheme, say MS. The first issue has been addressed by the community by simu- 
lating three different discretizations of the Dirac action - Wilson, SW clover, and staggered 
- which have discretization errors of 0(a), 0(a s (a)a), and 0(a 2 ) respectively. The second 
issue, reliability of the 1-loop perturbative matching relations, is being checked by using 
non-perturbative estimates. 

For Wilson and SW clover formulations, the internal consistency of the lattice calculations 
can be checked by calculating the quark masses two different ways. The first is based 
on methods of xPT, i.e. the calculated hadron masses are expressed as functions of the 
quark masses as in %PT. This method, based on hadron spectroscopy, is labeled HS for 
brevity. In the second method, labeled WI, quark masses are defined using the ward identity 
o '/LtV'75') = ( m i + m 2)iplbip- An example of such checks is shown in Fig. |2|. The solid lines 
are fits to the HS and WI estimates using the Wilson action and on the same statistical 
sample of configurations. The very close agreement of the extrapolated values is probably 
fortuitous since the linear extrapolation in a, shown in Fig. ^ neglects both the higher 
order discretization errors and the 0(a s (a) 2 ) errors in the 1-loop perturbative matching 
relations. The figure also shows preliminary results for the same WI data but now with 
non-perturbative Z's. The correction is large, note the large change in the slope in a, yet 
the extrapolated value is ^ 4 MeV. The final analysis using non-perturbative Z's will be 
available soon, and it is unlikely that the central value presented below will shift significantly. 

Lastly, in the quenched approximation, estimates of quark masses can and do depend on 
the hadronic states used to fix them. The estimates of fn given below were extracted using 
the pseudoscalars mesons, i.e. pions, with the scale a set by M p . Using either the nucleon or 
the A to fix m would give ~ 10% smaller estimates. Similarly, extracting m s using Mk gives 
estimates that are ~ 20% smaller that those using Mk* or as shown below. While these 
estimates of quenching errors are what one would expect naively, we really have to wait for 
sufficient unquenched data to quantify these more precisely. 

A summary of the quenched results in MeV at scale (i — 2 GeV, based on an analysis of 
the 1997 world data B2I, is 
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o CP-PACS R$ □ LANL U$ 
x CP-PACJg WI * LANL WI 
o Reweighted WI 
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Figure 2: Linear extrapolation of m versus a(M p ) for Wilson fermions using HS and WI 
methods. The WI data corrected by using non-perturbative estimates for the matching 
constants are also shown. 





Wilson 


TI Clover 


Staggered 


m(M n ) 
m s (M K ) 


4.1(1) 
107(2) 
139(11) 


3.8(1) 
99(3) 
117(8) 


3.5(1) 
91(2) 
109(5) . 



The difference in estimates between Wilson, tadpole improved (TI) clover, and staggered 
results could be due to the neglected higher order discretization errors and/or due to the dif- 
ference between non-perturbative and 1-loop estimates of Z's. This uncertainty is currently 



15%. Similarly, the 



i variation in m s with the state used to extract it, M K versus 



Mk* (or equivalently M$) could be due to the quenched approximation or again an artifact 
of keeping only the lowest order correction term in the extrapolations. To disentangle these 
discretization and quenching errors we again need precise unquenched data. 

Thus, for our best estimate of quenched results we average the data and use the spread as 
the error. To these, we add a second uncertainty of 10% as due to the determination of the 
scale 1/a (another estimate of quenching errors). The final results, in MS scheme evaluated 
at 2 GeV, are 1321 



3.8(4) (4) MeV 
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m s = 110(20)(11) MeV. 



(11) 



The important question is how do these estimates change on unquenching. The 1996 
analyses suggested that unquenching could lower the quark masses by ~ 20% [|4], |35|, how- 



ever, as discussed in I no longer feel confident making an assessment of the magnitude 
of the effect. The data does still indicate that the sign of the effect is negative, i.e. that 
unquenching lowers the masses. An estimate of the size requires more unquenched data. 

To end this section let me comment on a comparison of the quenched estimates with 
values extracted from sum-rules as there seems to be a general feeling that the two estimates 
are vastly different. In fact the recent analyses indicate that the quenched lattice results and 
the sum-rule estimates are actually consistent. A large part of the apparent difference is due 
to the use of different scales at which results are presented. Lattice QCD results are usually 
stated at = 2 GeV, while the sum-rules community uses \i = 1 GeV, and the running of 
the masses between these two scales is an « 30% effect in full QCD. This issue is important 
enough that I would like to briefly review the status of sum-rule estimates. 



6 Sum rule determinations of m and m 



s 



A summary of light quark masses from sum-rules is given in Table |l| Sum rule calculations 
proceed in one of two ways, (i) Using axial or vector current Ward identities one writes a 
relation between two 2-point correlation functions. One of these is evaluated perturbatively 
after using the operator product expansion, and the other by saturating with intermediate 
hadronic states |36|, [37]]. The quark masses are the constant of proportionality between these 
two correlation functions, (ii) Evaluating a given correlation function both by saturating 
with known hadronic states and by evaluating it perturbatively p8 |. The perturbative ex- 



pression depends on quark masses, and defines the renormalization scheme in which they are 
measured. The main sources of systematic errors arise from using (i) finite order calculation 
of the perturbative expressions, and (ii) the ansatz for the hadronic spectral function. Of 
these the most severe is the second as there does not exist enough experimental data to 
constrain the spectral function even for fi < 2 GeV. Since there are narrow resonances in 
this region, one cannot match the two expressions point by point in energy scale. The two 
common approaches are to match the moments integrated up to some sufficiently high scale 
(finite energy sum rules) or to match the Borel transforms. The hope then is that the result 
is independent of this scale or of the Borel parameter. 

Progress in sum-rules analyses has also been incremental as in LQCD. The perturbative 
expressions have now been calculated to 0{a z s ) |fH)| , and the value of J^q CD has settled at 
380MeV. A detailed analysis of the convergence of the perturbation expansion suggests 
that the error associated with the truncation at 0(a^) is ~ 10% for \i > 2 GeV . 



Improving the spectral function has proven to be much harder. For example, Colangelo 
et al. [[Hp have extended the analysis of m s in |J7|, |4T| by constructing the hadronic spectral 
function up to the first resonance (if* (1430)) from known Kir phase shift data. Similarly, 
Jamin 03 has used a different parametrization of the Omnes representation of the scalar 



form factor using the same phase shift data. In both cases the reanalysis lowers the estimate 
of the strange quark mass significantly. The new estimates, listed in Table [L| are consistent 
with the quenched estimates discussed in Section 
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reference 



m (MeV) 



m s (MeV) 



P| 1989 



|| 1995 
1 3|| 1995 
1 53 1995 
H 1996 



m 1997 



H 1997 
|43[ 1997 



|45| 1997 



4§ 1997 



47j 1997 



6.2(0.4) 
4.7(1.0) 
5.1(0.7) 



= 4.9(1.9) 
> 3.8-6 

> 3.4 
> 4.1-4.4 



= 138(8) 

= 144(21) 
= 137(23) 
= 148(15) 
= 91 - 116 
= 115(22) 

> 118 - 189 
> 88(9) 

> 104-116 



Table 1: Values and bounds on m and m s , in MS scheme at 2 GeV, from sumrule analyses. 



One can circumvent the uncertainties in the ansatz for the spectral function by deriving 
rigorous lower bounds using just the positivity of the spectral function |33], [15], |16|, fT7|| . Of 
these the most stringent are in [47| which rule out m < 3 and m s < 80 MeV for fi < 2.5 GeV. 
The bounds, however, have a significant dependence on the scale fi as evident by comparing 
the above values to those in the last row in Table [I], and the open question is how to fix fi, i.e. 
the upper limit of integration in the finite energy sum rules at which duality between PQCD 
and hadronic physics becomes valid? Unfortunately, this question cannot be answered ab 
initio. 



7 Implications for e'/e 

The Standard Model (SM) prediction of e'/e can be written as [18 

e'/e = A jc + [c 6 5 6 1/2 + c 8 5 8 3/2 ]M r j , (12) 

where M r = (158MeV/(m s + m^)) 2 and all quantities are to be evaluated at the scale 
m c = 1.3GeV. Eq. O highlights the dependence on the light quark masses and the bag 

1/2 3/2 

parameters B 6 and B s . For the other SM parameters that are needed in obtaining this 
expression we use the central values quoted by Buras et al. HSJ . Then, we get A = 1.29 x 10 -4 , 



c = —1.4, c 6 = 7.9, c 8 = —4.0. Thus, to a good approximation e'/e cx M r . 

Conventional analysis, with m s + = 158MeV and B^ 2 = B^ 2 = 1, gives e'/e ~ 
3.2 x 10~ 4 . The uncertainties in the remaining SM parameters used to determine A, cq,cq, 
and eg in Eq. [T2| are large enough that, in fact, any value between —1 x 10~ 4 and 16 x 10~ 4 is 
acceptable^pj . Current experimental estimates are 7.4(5.9) x 10 -4 from Fermilab E731 [|9 



and 23(7) x 10 4 from CERN NA31 [50]. So at present there is no resolution of the issue 
whether the CKM ansatz explains all observed CP violation. 

The new generation of experiments, Fermilab E832, CERN NA48, and DA$NE KLOE, 
will reduce the uncertainty to ~ 1 x 10~ 4 . First results from these experiments should be 
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available in the next couple of years. Thus, it is very important to tighten the theoretical 
prediction. 

As is clear from Eq. [12], both the values of quark masses and the interplay between B^ 2 
and B^ 2 will have a significant impact on e'/ e. The lower values of quark masses suggested 
by lattice QCD analyses would increase the estimate. The status of results for the various 
B-parameters relevant to the study of CP violation are discussed in the next section. 



8 B-parameters, Bk, Bq, B 7 , B 8 

Considerable effort has been devoted by the lattice community to calculate the various B- 
parameters needed in the standard model expressions describing CP violation. A summary 
of the results and the existing sources of uncertainties are as follows. 



8.1 B K 

The standard model expression for the parameter e, which characterizes the strength of the 
mixing of CP odd and even states in Kl and Ks, is of the form [[51] 

|e| ~ lm(V td V t :)BM H^r^iH) , (13) 

1V1\Y My/ 

where $ is a known function involving Inami-Lim functions and CKM elements. This rela- 
tion provides a crucial constraint in the effort to pin down the p and r\ parameters in the 



Wolfenstein parameterization of the CKM matrix [52]. The quantity Bk parameterizes the 
QCD corrections to the basic box diagram responsible for K° — K° mixing. This transition 
matrix element is what we calculate on the lattice. 

The calculation of Bk is one of the highlights of LQCD simulations. It was one of the 
first quantities for which theoretical estimates were made, using quenched chiral perturbation 
theory, of the lattice size dependence, dependence on quark masses, and on the effects of 
quenching ]53], [54], [55]]. Numerical data in the staggered formulation (which has the 



advantage of retaining a chiral symmetry which preserves the continuum like behavior of 
the matrix elements) is consistent with these estimates in both the sign and the magnitude. 
Since all these corrections have turned out to be small, results for Bk with staggered fermions 
have remained stable over the last five years. 

Three collaborations have pursued calculations of Bk using staggered fermions. The 
results are B K (NDR, 2 GeV) = 0.62(2)(2) by Kilcup, Gupta, and Sharpe |5|, 0.552(7) 



by Pekurovsky and Kilcup and 0.628(42) by the JLQCD collaboration 0. Of these, 
the results by the JLQCD collaboration are based on a far more extensive analysis. The 
quality of their data are precise enough to include both the leading 0(a 2 ) discretization 
corrections, and the 0{a 2 s ) corrections in the 1-loop matching factors. Their data, along 
with the extrapolation to a = limit including both factors, are shown in Fig. [|. I consider 
theirs the current best estimate of the quenched value. 

The two remaining uncertainties in the above estimate of Bk are quenching errors and 
SU(3) breaking effects (all current results have been obtained using degenerate quarks m u = 
r^d — m si with kaons composed of two quarks of mass ~ m s /2 instead of m s and m^). There 
exists only preliminary unquenched data K| which suggest that the effect of sea quarks 



13 



B K (NDR, 2GeV) vs. m p a 

q =1/a, 3-loop coupling, 5 points 




o non-invariant 
O invariant 



0.5 1 1 1 1 1 1 1 1 1 1 1 1 1 

-0.2 0.0 0.2 0.4 0.6 0.8 

m p a 

Figure 3: JLQCD data for Bk(NDR,2 GeV), and their extrapolation to the continuum 
limit keeping both the 0(a 2 ) discretization corrections, and the 0(a 2 ) in the 1-loop matching 
factors for two different discretizations of the weak operator. The extrapolation of the data 
without including the 0(a 2 ) corrections are shown by the dashed lines. 

is to increase the estimate by ~ 5%. Lastly, Sharpe has used x?T to estimate that the 
SU(3) breaking effects could also increase Bk by another 4 — 8% PSfl . Confirmation of these 
corrections requires precise unquenched data which is still some years away. 

In a number of phenomenological applications what one wants is the renormalization 
group invariant quantity Bk defined, at two-loops, as 

D D ( \ ( ( \\-7o/2/3 /i i a s{[l) r /3l70 — A)70i\ / 1/n 

B K = B K {ji) (a s (/i)) po (l + — — [ —2 ]) • (14) 

47T 2fj$ 

Unfortunately, to convert the quenched JLQCD number, 0.628(42), one has to face the issue 
of the choice of the value of a s and the number of flavors rif. It turns out that the two- loop 
evolution of Bk is such that one gets essentially the same number for the quenched theory, 
0.87(6) for rif = and a s (2GeV) = 0.192, and 0.84(6) for the physical case of rif = 3 and 
a s (M T ) = 0.354. One might interpret this near equality to imply that the quenching errors 
are small. Such an argument is based on the assumption that there exists a perturbative 
scale at which the full and quenched theories match. Since there is no a priori reason to 
believe that this is true, my preference is to double the difference and assign a second error of 
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0.06, an estimate also suggested by the preliminary unquenched data and the xPT analysis. 
With this caveat I arrive at the lattice prediction 



B K = 0.86(6)(6). 



(15) 



8.2 (7t + 7t°\O a \K + ) and its relation to B K 

The AS = 2 operator sL^dsL^d responsible for the K° — K° transition belongs to the same 
27 representation of SU(3) as the AS" = 1, AI = 3/2 operator C 4 = sL n d(uL^u — dL^d) + 
sL^uuL^d. At tree level in x?T, one gets the relation [^9], Q 



(^ U |Oa5= 2 |^°) 



2M| 



3i 



Ml 



Ml 



71 



~7c \O A \K + ) 



(16) 



So one way to calculate Bk is to measure (n + ii \O4\K + ) on the lattice and then use Eq. |16| . 
The motivation for doing this is that, for Wilson-like lattice actions, is only multiplica- 
tively renormalized due to CPS symmetry [61], whereas the AS = 2 operator mixes with all 
other chirality operators of dimension six. Using 1-loop values for these mixing coefficients 
has proven inadequate, though the recent implementation of non-perturbative estimates has 



made this situation much better 62 



The first lattice calculations of the AI = 3/2 part of the K + — > 7r + 7r° amplitude |63], 65 1 
gave roughly twice the experimental value, even though the Bk extracted from this "wrong" 
amplitude "agrees" with modern lattice estimates. The important question therefore is how 



reliable are the ;\TT relations between 



,7T 1 7T 



0,\K + )\ lattlce and (7r+7r°|e> 4 |if + >| 



between (7r + 7r |O 4 |^ 



physical' 



and 



Wattice 



and (K \O &s=2\K°) , i.e. does one or both fail? 



There are three sources of systematic errors in lattice calculations of \0^\K + ) that 

could explain the contradiction. The calculation is done in the quenched approximation, on 
finite size lattices, and with unphysical kinematics (in the lattice calculation the final state 
pions are degenerate with the kaon since m u = irtd = m s , and are at rest). Of the three 
possibilities, the last is the most serious as it gives a factor of two even at the tree-level 



Ml 



2m2r {h + AOa\K- 



Ml 



\ physical 



(7l + 7T°\0 4 \K A 



I unphysical 



(17) 



This tree-level correction was taken onto account in [^3| |65|. The source of the remaining 



discrepancy by a factor of two was anticipated by Bernard as due to the failure of the tree- level 
expression 0]. Recently, Golterman and Leung []59f have calculated the 1-loop corrections 
to Eq. 17. This calculation involves a number of unknown 0(p A ) chiral coefficients of the 



weak interactions and thus has a number of caveats. However, under reasonable assumptions 
about the value of these 0(p 4 ) constants, the corrections due to finite volume, quenching, and 
unphysical kinematics all go in the right direction, and the total 1-loop correction can modify 
Eq. [L7| by roughly a factor of two. On the other hand the modification to the connection 
with Bk at the physical point is small. 

The JLQCD Collaboration |66| has recently updated the calculations in |65|]. By 
improving the statistical errors they are able to validate the trends predicted by 1-loop 
expressions. Thus one has a plausible resolution of the problem. I say plausible because the 
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calculation involves a number of unknown chiral couplings in both the full and quenched 
theory and also because of the size of the 1-loop correction. The conclusive statement is the 
failure of xPT for the relation Eq. [17| . 

The other relevant question is what bearing does the analysis of Golterman and Leung 



59| have on the calculations of other S-parameters. In the calculation of Bk, based on 
measuring the (K°\Oas=2\K°) transition matrix element as discussed in section |8.1| , ;\PT 
has been used only to understand effects of finite volume and chiral logs. These are found to 
be small and the data show the predicted behavior. Based on this success of xPT we estimate 
that the two remaining errors - quenching and the use of degenerate quarks - are each « 5% 

3/2 

as suggested by 1-loop ;\PT. On the other hand, in present calculations of B$, B 7 , and 
B^ 2 xPT is used in an essential way, i.e. to relate (7r + ir \O\K + } to (tv + \0\K + } . Second, 
the calculations are done for unphysical kinematics, i.e. the final state pion is degenerate 
with the kaon. It would be interesting to know the size of the one-loop corrections to these 
relations. 



8.3 Sy 72 and bI' 2 

Assuming the reduction of {-R + ^\Oji =z ' 2 \K + ) to (tt+| Cfg =3/2 |ir + > using tree-level xPT is 
reliable, the lattice calculations of B^ 2 and B^ 2 are as straightforward as those for Bk- 
There are three "modern" quenched estimates of B^ 2 and B% 2 which supercede all previous 
reported values. These, in the NDR-MS scheme at 2 GeV, are 



Fermion type 


Matching Z 




l/a GeV 


Bp 2 


„3/2 


(A) Staggered [55| 


1-loop 


6.0, 6.2 


a -> 


0.62(3)(6) 


0.77(4) (4) 


(B) Wilson |7| 


1-loop 


6.0 


2.3 


0.58(2) (7) 


0.81(3)(3) 


(C) Tree-level Clover 62 




1-loop 


6.0 


2.0 


0.58(2) 


0.83(2) 


(D) Tree-level Clover [52 




Non-pert. 


6.0 


2.0 


0.72(5) 


1.03(3) . 



where I have also given the type of lattice action used, the /?'s at which the calculation was 
done, the lattice scale l/a at which the results were extracted, and how the 1-loop matching 
coefficients were determined. 

The difference between (C) and (D) is the use of perturbative versus non-perturbative 
Z's. Thus (D) is the more reliable of the APE numbers. The agreement between (B) and 
(C), in spite of the difference in the action, is a check that the numerics are stable. All three 
of these results suffer from the fact that these calculations were done at (3 = 6.0 (l/a « 2 
GeV) and there does not yet exist data at other f3 needed to do the extrapolation to the 
continuum limit. 

The result (A) does incorporate an extrapolation to a=0, but with only two beta values. 
For example, B^ 2 = 1.24(1) and 1.03(2) at (3 = 6.0 and 6.2 respectively. Due to the large 
slope in a 2 , such an extrapolation based on two points should be considered preliminary. 
Lastly, one needs to demonstrate that corrections to 1-loop Z's are under control. 
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8.4 B 6 



The recent work of Pekurovsky and Kilcup provides the best lattice estimate for Bq. 
Their results Bq = 0.67(4) (5) for quenched and 0.76(3) (5) for two flavors have the following 
systematics that are not under control. The calculation is done for degenerate quarks, 
m u = md — m s and uses the lowest order to relate (irir\0\K + ) to (tt\0\K + ). There 
is no reason to believe that higher order corrections may not be as or more significant as 
discussed above for C 4 . The second issue is that the 1-loop perturbative corrections in the 
matching coefficients are large. Lastly, there is no estimate for the discretization errors as 
the calculation has been done at only one value of (3. Thus, at this point there is no solid 
prediction from the lattice. 



9 Conclusions and Acknowledgements 

In view of the new generation of ongoing experiments to measure e'/e with the proposed 
accuracy of 1 x 10 -4 , it is very important to firm up the theoretical prediction. The standard 
model estimate depends very sensitively on the sum m s + m d and on the interplay between 
the strong and electromagnetic penguin operators, i.e. B X J 2 and B^ 2 . Quenched lattice 
results for (m s + m^)(2 GeV) are settling down at 115(25) MeV, and preliminary evidence is 
that unquenching further lowers these estimates. The calculations of B X J 2 and B^ 2 are less 
advanced. Hopefully we can provide reliable quenched estimates for these parameters in the 
next year or so. Thereafter, we shall start to chip away at realistic full QCD simulations. 

Acknowledgement s 

I would like to thank the organizers of CPMASS97 and ORBIS SCIENTIAE'97 for 
very stimulating conferences, and my collaborators Tanmoy Bhattacharya, Greg Kilcup, 
and S. Sharpe for many discussions. I also thank the DOE and the Advanced Computing 
Laboratory for support of our work. 



References 



M. Creutz, Quarks, Gluons, and Lattices, Cambridge University Press, 1983 
Quantum Fields on the Computer, Ed. M. Creutz, World Scientific, 1992. 
Quantum Fields on a Lattice, I. Montvay and G. Minister, Cambridge Monographs on 
H. J. Rothe, Quantum Gauge Theories: An Introduction, World Scientific, 1997. 
H.B.Nielsen and M. Ninomiya, [Nucl. Phys. B185 (1981T20 . 



Y. Shamir, |Nucl. Phys. (Proc. Suppl.) B53 (19977661] ; [INucl. Phys. (Proc. Suppl.) B47 
(1996) 212j . 



R. Narayanan, [hep- lat/ 9 70 70 3 5 



M. Testa, |hep-lat/97070TH . 

M. Alford, T. Klassen, and P. Lepage, |Nucl. Phys. B496 (1997) 377fc and [Eep^ 
lat/9712005 



17 



M. Liischer, et al, [Nucl. Phys. B478 ( 1996T365 
P. Hasenfratz, |hep-Iat/9709TTU 



K. G. Wilson, phys. Rev. DIP (1974) 2445| ; in New Phenomena in Subnuclear Physics 
(Erice 1975), Ed. A. Zichichi, Plenum New York, 1977. 



B. Sheikholesalami and R. Wohlert, |Nucl. Phys. B259 (1985T572 



J. Kogut and L. Susskind, |Phys. Rev. Dll (1975) 395| ; T. Banks, J. Kogut, and 
L. Susskind, |Phys. Rev. D13 (1996) 1043]; L. Susskind, |Phys. Rev. Dl6 (1977) 3031. 



M. Liischer, Selected Topics In Lattice Field Theory, Lectures at the 1988 Les Houches 
summer school, North Holland, 1990. 

K. Symanzik, in Mathematical Problems in Theoretical Physics, Springer Lecture Notes 
in Physics, bf 153 (1982) 47; M. Liischer and P. Weisz, it Commun. Math. Phys. 97 
1985 59; G. Heatlie, et al, |Nucl. Phys. B352 (1991)"266 . 



H. Wittig, |hep-lat/97100T3| ; M. Peardon, |hep-lat/9710029 
P. Lepage and P. Mackenzie, Phys. Rev. D48 (1993) 225C| . 

P. Lepage and B. Thacker, |Phys. Rev. D43 (1991)~T96| ; G. Lepage, L. Magnea, C. 
Nakhleh, U. Magnea and K. Hornbostel, Phys. Rev. D46 (1992) 405% 

A. El-Khadra, A. Kronfeld, and P. Mackenzie, phys. Rev. D55 (1997) 3933 . 



M. Neubert, phys. Rev. 245 (1994) 259| ; G. Martinelli and C. Sachrajda, |Nucl. Phys. 
B478 (1996) BSD 



J. Mandula and M. Ogilvie, [hep- lat/ 9 70 3020] . 



G.C. Rossi, |lNucl. Phys. (Proc. Suppl.) B53 (197Tj~3] M. Liischer, et al, [Nucl. Phys 



B478 (1996)~365l G. Martinelli, et al, [Nucl. Phys. B445 (1995 [81 



S. Sharpe, |Nucl. Phys. (Proc. Suppl.) B17 (1990) 146[ Phys. Rev. D41 (1990) 1990 



Phys. Rev. D46 (1992) 31461 , |Nucl. Phys. (Proc. Suppl.) B30 (1993HT3 . 
C. Bernard and M. Golterman, phys. Rev. D46 (1992)~g53]; [Nucl. Phys. (Proc. SuppT 



B30 (1993) ~2T7| ; [Nucl. Phys. (Proc. Suppl.) B37 (1994) 334! ; M. Golterman, [hep 



lat/9405002| and |hep-lat/9411005 



R. Gupta, |lNucl. Phys. (Proc. Suppl.) B42 (1995)~85| . 
S. Sharpe, [Nucl. Phys. (Proc. Suppl.) B53 (19977181 



J. Gasser et al, phys. Rep. C87 (1982) 77 . 
H. Leutwyler, [INucl. Phys. B337 (1990) T08 



H. Leutwyler, [hep-ph/ 9609467 . 

D. Kaplan and A. Manohar, phys. Rev. Lett. 56 (1986) 2004 . 
T. Banks, Y. Nir, N. Sieberg, pep-ph/9403203 . 
T. Bhattacharya, R. Gupta, |hep- lat/9710095| . 

T. Bhattacharya, R. Gupta, G. Kilcup, and S. Sharpe, phys. Rev. D53 (1996) 6486 . 
R. Gupta and T. Bhattacharya et al, Phys. Rev. D55 (1997) 7203 



B. Gough et al, |Phys. Rev. Lett. 79 (1997P522 



18 



J. Bijnens, J. Prades, and E. de Rafael, Phys. Lett. 348B (1995l~226 
M. Jamin and M. Mlinz, |Z. Phys. C66 (1995) 633 . 



S. Narison, |Phys. Lett. 358B (1995)TT3 



S. Narison, |Phys. Lett. 216B (1989) 335 



K.G. Chetyrkin, D. Pirjol, and K. Schilcher, Phys. Lett. 404B (1997)1337 . 
P. Colangelo, F. De Fazio, G. Nardulli, and N. Paver, Phys. Lett. 408B (1997T340 



M. Jamin, |hep-ph/9709484 . 
J. Prades, |hep-ph/9708395| . 

T. Bhattacharya, R. Gupta, K. Maltman, [hep-ph / 9703455 . 
F.J. Yndurain, |hep-ph/9708"300 



H.G. Dosch and S. Narison, [hep-ph/9709215 



L. Lellouch, E. de Rafael, J. Taron, [Phys. Lett. 414B (1997)^95 



A. Buras, M. Jamin, and M. E. Lautenbacher, |Phys. Lett. 389B (1996) 749| ; [hej> 
ph/9608365 



L. K. Gibbons, et ai, |Phys. Rev. Lett. 70 (1993) 1203 



G. D. Barr, et ai, |Phys. Lett. 317B (1993)~233 



G. Buchalla, A. Buras, and M. E. Lautenbacher, |Rev. Mod. Phys. 68 (1996) 1125 
A. Buras, |hep-ph/9711217| . 



G. Kilcup, S. Sharpe, R. Gupta, and A. Patel, |Phys. Rev. Lett. 64 (1990) 25 



S. Sharpe, [Nucl. Phys. (Proc. Suppl.) B34 (19947403 



G. Kilcup, R. Gupta, and S. Sharpe, |Phys. Rev. D57 (1998) 1654] |hep-lat/ 9707001 . 
D. Pekurovsky and G. Kilcup, |hep-lat/9709T46 



S. Aoki, et ai, JLQCD collaboration, [hep-lat/9710073 . 

G. Kilcup, D. Pekurovsky, and L. Venkatraman, |Nucl. Phys. (Proc. Suppl.) B53 (1997) 
J4§ 



M. Golterman and K. C. Leung, |Phys. Rev. D56 (1997) 2950| . 
J. Donoghue, E. Golowich, and B. Holstein, [Phys. Lett. 119B (1982")"4T2 . 
C. Bernard, et ai, [Phys. Rev. D32 (1985) 2~343| . 
L. Conti, A. Donini, V. Gimenez, G. Martinelli, M . Talevi, A. Vladikas, [hep-lat/9711053 



C. Bernard and A. Soni, [Nucl. Phys. (Proc. Suppl.) B9 (1990TT55| ; [INucl. Phys. (Proc 



Suppl.) B17 (1990) 495 



C. Bernard, in From Actions to Answers, Proceedings of the 1989 TASI School, edited 
by T. DeGrand and D. Toussaint, (World Scientific, 1990). 

M.B.Gavela, et ai, [Phys. Lett. 211B (1988) gg ; |Nucl. Phys. B306 (1988) 677 . 

S. Aoki, et a/., JLQCD Collaboration, |hep-lat/971IM6" , 



R. Gupta, T. Bhattacharya, and S. Sharpe, [Phys. Rev. D55 (1997) 403T) 
G. Kilcup, |lNucl. Phys. (Proc. Suppl.) B20 (1991J4T7 . 



19 



